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A search for lepton flavor violating decays of the r lepton to a lighter mass lepton and a pseu- 
doscalar meson has been performed using 339 fb -1 of e + e~ annihilation data collected at a center-of- 
mass energy near 10.58 GeV by the BABAR detector at the SLAC PEP-II storage ring. No evidence 
of signal has been found, and upper limits on the branching fractions are set at 10 -7 level. 

PACS numbers: 13.35.Dx, 14.60.Fg, 11.30.Hv 



The recent discovery of large neutrino mixing [lj sug- 
gests that lepton flavor violation (LFV) occurs. Charged 
LFV decays have not yet been observed, although they 
have long been identified as unambiguous signature of 
new physics. Ncutrinolcss decays like r ± — > £ ± P°, where 
I = e,fi and P Q = 7r°, 77,77', are likely candidates for 
LFV 0,13, which could be induced by potentially large 
mixing between the supersymmetric partners of the lep- 
tons and is further enhanced by color factors associ- 
ated with these semi-leptonic decays. Some models with 
heavy Dirac neutrinos 0, @|, two Higgs doublet mod- 
els, R-parity violating supersymmetric models, and fla- 
vor changing Z' models with non-universal couplings JfJ 
allow for observable parameter space of new physics [7|, 
while respecting the existing experimental bounds (3 • 

The results presented here use an integrated luminosity 
L= 339 fb -1 collected at a center-of-mass (CM) energy, 
y/s, near 10.58 GeV by the BABAR detector at the SLAC 
PEP-II e + e~ asymmetric-energy storage ring. Details of 
the BABAR detector are described elsewhere [9| . 

The signature of the signal process is the presence of an 
£P° pair having an invariant mass consistent with m T 
1 .777 GeV/c 2 [l(| and a total energy equal to y/s/2 in the 
CM frame, along with other particles in e + e~ — > t + t~ 
events having properties consistent with a r lepton de- 
cay. Two neutral decay modes (ir° — > 77 and 77 — > 77) 
and three charged decay modes [77 — > 7r + 7r~7r° (tt° — > 77), 
77' — > 7r + 7r~?7 (77 — > 77), and 77' — > p°7 ] are reconstructed. 
Signal events are simulated with the KK2F [llj Monte 
Carlo (MC) program, where the r — > ^ ± P° decays ac- 
cording to two body phase space, while the other r decays 
according to measured branching fractions [l2| simulated 
with TAUOLA [13 . /j + \x~ and r + r~ background processes 
are generated using KK2F and TAUOLA, and qq processes 
are generated using EVTGEN [lj] and JETSET [15]. Ra- 
diative corrections are simulated using PHOTOS |la| . The 
detector response is simulated with GEANT4 17]. The 
MC events are used for the optimization and systematic 
studies of the signal efficiencies, and for determination of 
the background shapes. Estimates of the rates for the 
backgrounds are derived directly from the data. 

Events with two or four well reconstructed tracks and 
zero total charge are selected. Tracks are rejected if 
they are consistent with coming from photon conversions. 



An event is divided into two hemispheres ("signal"- and 
"tag"- sides) in the CM frame by a plane perpendicu- 
lar to the thrust axis [13, calculated using all observed 
particles. 

The signal-side hemisphere is required to contain one 
or three tracks and two photon candidates with energy 
E 1 > 50MeV for the tt° -> 77, 77 -> tt+tt^tt (tt° -> 77) 
and 77' — > 7r + 7r~?7 (77 — » 77) channels, and E 1 > 100 MeV 
for the 77 — * 77 channel. For the 77' — * p°7 channel, 
the single photon candidate is required to have E y > 
100 MeV. Events with additional photon candidates in 
the signal hemisphere with E 1 > 100 MeV are rejected. 

The P° candidates are reconstructed in the following 
mass windows: m(7r° — * 77) e [0.115, 0.150] GeV/c 2 , 
771(77 -> 77) S [0.515, 0.565] GeV/c 2 , 777(77 tt+tt-tt ) 
€ [0.537, 0.558] GeV/c 2 , 777(77' -> ^ + tt'v) G [0-950, 
0.965] GeV/c 2 , 777(77' -> P°l) € [0.940, 0.970] GeV/c 2 , 
and 777(p° -> tt+tt-) e [0.600, 0.900] GeV/c 2 . To reduce 
combinatorial backgrounds, the momentum of P° is re- 
quired to satisfy: p^o > 0.5GeV/c for — > e ± -K Q , p^o > 
1.5GeV/c for t± -> ^tt , p v > l.OGeV/c for r± e ± r t 
(77 -> 77), p n > 1.4 GeV/c for -> p^r] (77 77), p v > 
1.2GeV/c for —> /i ± ?7 (77 — * ir + TT~ir°, n° — > 77) and 
p,,' > 1.3GeV/c for t ± — * £^77' (77' — > p°7) decays. 

The track unassociated with any of the P° daughters is 
required to have a momentum > 0.5 GeV/c and is identi- 
fied as an electron or muon, but not as a kaon, using stan- 
dard BABAR. particle identification techniques [l9( . In the 
case of a charged P° decay, this criteria is applied on the 
track that combines with the one having opposite-sign 
charge and the photon candidate (s) to give an invariant 
mass farthest from the nominal P° mass [13 • This pro- 
vides the correct pairing for > 99.7% of selected signal 
MC events after particle- identification requirements. 

The origin of the photon(s) is assigned to the point of 
closest approach of the lepton track to the e + e _ collision 
axis for neutral P° decays, or to the common vertex in the 
signal-side hemisphere for the charged P° decays. The 
P° momentum is kinematically fitted with its respective 
mass constraints, and combined with the lepton track to 
form the signal r candidate. An event is accepted based 
upon the closeness of the signal r candidate to 777 T . 

Signal decays are identified by two kinematic vari- 
ables: the beam-energy constrained t mass 777ec and 
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AE = Ee+Epa where Eg and Epo are the respec- 

tive energies in the CM frame. These two variables are 
independent apart from small correlations arising from 
initial and final state radiation. For signal events, the re- 
constructed peak positions of the mEC distribution agree 
very well with m T , while those of AE vary between —5 
to — 23MeV. The shift from zero in the AE peak comes 
from mis-calibration of the measured photon energy. The 
resolutions of the ttiec and AE distributions for the sig- 
nal events are presented in Table [J Events in the data 
within a ±3cr rectangular box centered around the signal 
MC peak positions are excluded until all optimization 
and systematic studies of the selection criteria have been 
completed. The selection i s op timized to yield the small- 
est expected upper limit [20| in a background-only hy- 
pothesis for observing events inside the ±2tr rectangular 
signal box around the signal MC peak positions shown 
in Fig. |U 

The dominant backgrounds are from r — * evVj or 
t — ► pv decays in t + t~ events, with additional contri- 
butions from Bhabha, di-muon and qq processes. The 
backgrounds are higher for searches with muons, due to 
misidentification of a it track as a fi candidate. Another 
source of background is the mis-reconstruction of 77 and 
7T° candidates. 

Non-r backgrounds with radiation along the beam di- 
rections are suppressed by requiring the polar angle of the 
missing momentum to lie between —0.76 and 0.92. The 
total CM momentum of all tracks and photon candidates 
on the tag-side is required to be less than 4.75 GeV/c. 

A tag-side hemisphere containing a single track is clas- 
sified as e-tag, /J-tag or /i-tag if the track is exclusively 
identified as an electron, muon, or neither, respectively. 
For these tags, the total neutral CM energy in the hemi- 
sphere SE^ M is required to be less than 0.2 GeV, and 
the invariant mass mtag, calculated using all observed 
charged and neutral particles, to be less than 0.4GeV/c 2 . 
For — > e ± P° channels, the data events in e-tag are 
used as a control sample to estimate the Bhabha back- 
ground, and are not included in the final selection. If 
the track is neither an electron nor a muon, SE^ M > 
0.2 GeV and m tag € [0.6, 1.3] GeV/c 2 , the event is classi- 
fied as a p-tag. For searches of neutral P° decay modes, 
events with three tracks in the tag-side with m ta g € [0.9, 
1.6] GeV/c 2 are also allowed. 

Taking the direction of the tag-side r to be opposite 
the signal candidate, all tracks and photon candidates 
in the tag-side hemisphere are used to calculate the in- 
variant mass squared of the tag-side missing momentum 
(m 2 ). To reduce non-r backgrounds for — * e^ir , 
— > e ± ?7 (77 — * 77) searches, (m 2 /1.8 GeV 2 /c 4 ) — 
ln(2 xp^ iss / v / s)/2.0 is required to be less than unity [2 if . 
where p^iss i s the component of the missing momentum 
transverse to the collision axis. For the other searches, 
— ln(2 x Pmiss/v 7 ^) i s required to be less than 2.5, ex- 
cept for t~^~ — » e^i] (77 — > 7r + 7T~7r°) and — > fi Tjf 



(rj' — > Tr + iT~ri) searches, where very few events are ex- 
pected. 

To focus on selected signal-like events, a Grand Side 
Band (GSB) is defined in the toec vs. AE plane as: ttiec 
€ [1.5,2.0] GeV/c 2 and AE e [-0.8,0.4] GeV. With 
electrons as the lepton track, 22, 18, 4, 1 and 30 events 
survive in the GSB for tt° — > 77, 77 — * 77, 77 — * 7r + 7r~7r°, 
77' — * 7r + 7r~?7 and 77' — > p°j channels, and 311, 69, 24, 
24 and 285 events survive for the corresponding channels 
with muons as the lepton track, as shown by dots in 
Fig. [TJ Also shown are the shaded regions containing 
68% of the selected signal MC events inside the GSB. 

The number of expected background events in the 
signal box is extracted from an unbinned maximum 
likelihood fit to the distributions of toec and AE in 
data inside the non-blinded parts of the GSB, using 
two-dimensional probability density functions (PDF) for 
e + e~, t + t~ and qq backgrounds. The kernel of 

the PDFs are estimated [22j from the data control sam- 
ples for e + e~, and respective MC events for the others. 

The dominant contribution to the uncertainty in back- 
ground estimation arises from the statistical precision on 
the selected data sample inside the non-blinded parts of 
the GSB, or from the variation of background compo- 
nents within ±lcr from their fitted values. The observed 
and expected events from the fit to the data inside the 
±3(7 to ztllcr annular boxes and the signal boxes are 
shown in Table [II which confirm good modeling of the 
backgrounds in data and show no evidence of signal. 

The largest systematic uncertainties in the signal re- 
construction efficiency are due to the signal track mo- 
mentum and the photon energy scale and resolution, es- 
timated by varying the peak position and resolution of 
the ttiec and AE distributions. The errors associated 
with the modeling of each selection variable are estimated 
from the relative change in signal efficiency when vary- 
ing the selection criteria by the difference between the 
data and MC events in the mean of that variable. Other 
sources of systematic uncertainties include those arising 
from trigger inefficiencies, tracking and neutral energy re- 
construction efficiencies, the signal lepton identification, 
beam-energy scale and spread, luminosity estimation and 
e+e- t+t- cross-section (a TT = 0.89 ± 0.02 nb) [H]. 
About 2.4 million MC events are used per channel, re- 
sulting in a negligible systematic uncertainty due to MC 
statistics. Although the signal MC events have been 
modeled using a two body phase space model, the re- 
sults obtained in this analysis are insensitive to this as- 
sumption as demonstrated by considering the two ex- 
treme cases of a V — A and a V + A form of interaction. 

The upper limits for — > ^ ± P° decays are calculated 
using Bfj° L = NlP L /(2Lcr TT Be) 1 where N$° L is the 90% 
confidence level (C.L.) upper limit on the number of sig- 
nal events inside the signal box, B and e are the branch- 
ing fraction [T3 ] and reconstruction efficiency of the signal 
decay mode under consideration. To obtain a combined 
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FIG. 1: Selected data (dots) and 68% of signal MC events (shaded region) inside the GSB region, and the ±2cr signal box. 
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TABLE I: The ttiec and AE resolutions for the signal MC events, the number of observed (obs.) and expected (exp.) events 
inside ±3cr to ±llcr boxes and ±2er box, the branching fractions (B), the efficiencies (e), and the 90% C.L. upper limits (UL). 



Decay modes 


o-(rriEc) 


a(AE) 


±3(j to ±llcr box 


±2cr box 


B 




UL (xlO" 7 ) 


MeV/c^ 


MeV 


obs. 


exp. 


obs. 


exp. 


(%) 


(%) 


obs. 


exp. 


r± -> e ± 7T° (tt -> 77) 


9.1 


46.4 


4 


5.37±1.14 





0.17±0.04 


98.80±0.03 


2.83±0.25 


1.3 


1.4 


r* -> /^tt (t" -» 77) 


9.0 


46.4 


43 


40.68±4.32 


1 


1.33±0.15 


98.80±0.03 


4.75±0.37 


1.1 


1.1 


T"^~ — » 6*77 {V ~ > 77) 


8.5 


42.6 


4 


4.99±1.18 





0.20±0.05 


39.38±0.26 


3.59±0.24 


2.5 


2.8 


T"*~ — > 6*77 (77 — » 7T + 7r _ 7r u ) 


5.9 


31.4 





0.64±0.32 





0.02±0.01 


22.43±0.40 


3.17±0.32 


5.4 


5.5 


— > 6*77 













0.22±0.05 


Be = 2.12±0.20 (%) 


1.6 


1.9 


r* -> ^77 (77 -> 77) 


8.3 


40.8 


20 


17.36±2.12 


1 


0.67±0.08 


39.38±0.26 


7.03±0.53 


1.9 


1.6 


— > /^ ± 77 (77 — > 7r + 7r^7r") 


5.6 


31.0 


3 


2.01±0.41 





0.08±0.02 


22.43±0.40 


3.67±0.32 


4.5 


4.8 


r* ^77 










1 


0.75±0.08 


Be = 3.59±0.41 (%) 


1.5 


1.3 


— > 6*77' (77' — » n + n~r]) 


5.9 


31.0 





0.14±0.14 





0.01±0.01 


17.52±0.56 


3.75±0.27 


5.8 


5.9 


± -kill \ 
r -> e 77 (77 -> p 7) 


4.4 


24.3 


2 


2.97±0.54 





0.11±0.03 


29.40±0.90 


2.98±0.28 


4.2 


4.5 















0.12±0.03 


Be = 1.53±0.16 (%) 


2.4 


2.6 


r* — > 77' (77' — > n + ir~r]) 


5.6 


29.1 


1 


2.42±0.47 





0.07±0.02 


17.52±0.56 


5.87±0.46 


3.6 


3.8 


± ± / / / \ 
T -»• /J 77 (77 -> p 7) 


4.1 


23.1 


13 


11.06±0.65 





0.42±0.03 


29.40±0.90 


3.90±0.46 


2.7 


3.7 


± ± / 
r -» /j 77 













0.49±0.04 


Be = 2.18±0.26 (%) 


1.4 


2.0 



upper limit with 77 and 77' decays, the observed and ex- 
pected background events and the signal efficiencies are 
added using Be — (B\ x E\ ± B2 x e 2 ), where Bi, B2 are 
the respective branching fractions and e\ and £2 are the 
corresponding efficiencies. This combination takes into 
account correlated uncertainties from the track and neu- 
tral cluster reconstruction efficiency and the signal lepton 
identification. The observed and the expected upper lim- 
its at 90% C.L. are presented in Table U including all con- 
tributions from systematic uncertainties 0, Hfl . These 
limits present up to a factor of four improvement over 
the previously published results [8|], except for r ± — * ^rj 



search, where the limit is similar. 

Mixing between left-handed smuons and staus allows 
one to translate the — > limit to an exclusion 

plot in the tan/3 vs. tua plane where tan/3 is the 
ratio of the vacuum expectation values of the two Higgs 
doublets and tua is the mass of the pseudoscalar Higgs. 
The excluded regions at 95% C.L. from this — * /i ± ?7 
search (< 1.9 x 10~ 7 ) with right-handed neutrino mass 
= 10 14 GeV/c 2 introduced via the seesaw mechanism are 
shown in Fig. [JH This result is competitive with those 
obtained from the direct searches for Higgs — > bb, t + t~ 
decays by CDF [26| and DO 27 1, and complementary 
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FIG. 2: Excluded regions in tan/3 vs. tua plane (see text). 
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to the region excluded by the LEP experiments with a 



top quark mass of 174.3 GeV/c 2 [28|, for two common 
scenarios of stop- mixing benchmark models [2{|: m™ ax 
and no-mixing obtained with the Higgs mass parameter 
/i = — 200GeV/c 2 shown by darker and lighter shaded 
regions respectively. 
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